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Abstract 

We have calculated on the computer the sum Bm of reciprocals of all 
47 known Mersenne primes with the accuracy of over 12000000 decimal 
digits. Next we developed Bm into the continued fraction and calcu- 
lated geometrical means of the partial denominators of the continued 
fraction expansion of Bm- We get values converging to the Khinchin's 
constant. Next we calculated the n-th square roots of the denominators 
of the n-th convergents of these continued fractions obtaining values ap- 
proaching the Khinchin-Levy constant. These two results suggests that 
the sum of reciprocals of all Mersenne primes is irrational, supporting 
the common believe that there is an infinity of the Mersenne primes. For 
comparison we have done the same procedures with slightly modified 
set of 47 numbers obtaining quite different results. Next we investigated 
the continued fraction whose partial quotients are Mersenne primes and 
we argue that it should be transcendental. 
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1 Introduction. 

The Mersenne primes M n are primes of the form 2 P — 1 where p must be a prime, see e.g. 
Sect. 2. VII]. The set of Mersenne primes starts with Mi = 2 2 - 1, M 2 = 2 3 - 1, M 3 = 2 5 - 1 
and only 47 primes of this form are currently known, see Great Internet Mersenne Prime 
Search (GIMPS) at www.mersenne.org. The largest known Mersenne prime has the value 
Mn = 2 43112609 - 1 = 3.164 7 2 6 9 33 . . . x io 12978188 . In general the largest known primes are 
the Mersenne primes, as the Lucas-Lehmer primality test applicable only to numbers of the 
form 2 P — 1 needs a multiple of p steps, thus the complexity of checking primality of M n is 
0(log(M n ))- Let us remark that algorithm of Agrawal, Kayal and Saxena (AKS) for arbitrary 
prime p works in about (9(log 7 ' 5 (p)) steps and modification by Lenstra and Pomerance has 
complexity C(log 6 (p)). 

There is no proof of the infinitude of M n , but a common belief is that as there are pre- 
sumedly infinitely many even perfect numbers thus there is also an infinity of Mersenne primes. 
S. S. Wagstaff Jr. in [2H] (see also J2Tj §3.5]) gave heuristic arguments, that M n grow doubly 
exponentially: 

log 2 log 2 M n ~ ne" 7 , (1) 

where 7 = 0.57721566... is the Euler-Mascheroni constant. In the Fig. [I] we compare the 
Wagstaff conjecture with all 47 presently known Mersenne primes M n . Of these 47 known 
M n = 2 P — 1 there are 27 with p mod 4 = 1 and 19 with p mod 4 = 3. It is in opposite to the 
set of all primes where the phenomenon of Chebyshev bias is known: for initial primes there 
are more primes p = 3 (mod 4) than p = 1 (mod 4), [T9] . 

In this paper we are going to exploit two facts about the continued fractions to support 
the conjecture on the infinitude of Mersenne primes: the existence of the Khinchin constant 
and Khinchin-Levy constant. We calculate the sum of reciprocals of the Mersenne primes 
Bm = SnV-Mnj if there is infinity of Mersenne primes then this number Bm should be 
irrational (at least, because it is probably even transcendental, as it is difficult to imagine the 
polynomial with some mysterious integer coefficients whose one of roots should be Bm)- 

There exists a method based on the continued fraction expansion which allows to detect 
whether a given number r can be irrational or not. Let 

r = [oo(r); ai(r), a 2 (r), a 3 (r), . . .] = a (r) H (2) 

ai(r) H 



1 

a2{r) H 

a 3 (r) + '•• 

be the continued fraction expansion of the real number r, where ao(r) is an integer and all a^{r) 
with k > 1 are positive integers. The quantities afc(r) are called partial quotients or the partial 
denominators. Khinchin has proved [13], see also [20], that 

lim (aUr) . . . aJr)) » = TT I 1 + — } =K^ 2.685452001 (3) 

n ^°° ^Jil m(m + 2) J v ; 

is a constant for almost all real r [HI §1.8] (the term clq is skipped in ([3])). The exceptions are of 
the Lebesgue measure zero and include rational numbers, quadratic irrationals and some irra- 
tional numbers too, like for example the Euler constant e = lim„_ >00 (l + = 2.7182818285 . . . 
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for which the n-th geometrical mean tends to infinity like y/n, see pm §14.3 (p. 160)]. The 
constant K is called the Khinchin constant. If the sequence 

K(r; n) = (a 1 (r)o 2 (r) . . . a n (r)) » (4) 

for a given number r tend to if for n — > oo we can regard it as an indication that r is irrational — 
all rational numbers have finite number of partial quotients in the continued fraction expansion 
and hence starting with some no for all n > no will be a n = 0. It seems to be possible to 
construct a sequence of rational numbers such that the geometrical means of partial quotients 
of their continued fraction will tend to the Khinchin constant. 

The Khinchin — Levy's constant arises in the following way: Let the rational P n (r) / 'Q n (r) 
be the n-th partial convergent of the continued fraction of r: 

P n (r) 

n "/ x = [oo(r); at(r), a 2 (r), a 3 (r), . . . , a n (r)]. (5) 

In 1935 Khinchin [12] has proved that for almost all real numbers r the denominators of the 
finite continued fraction approximations fulfill: 

lim (Q n (r)) 1/n = lim (L(r;n) = L (6) 

and in 1936 Paul Levy [H] found an explicit expression for this constant L: 

lim K/Q n {r) = e" 2/121og(2) = L — 3.27582291872 ... (7) 

L is called the Khinchin — Levy's constant [SJ §1.8]. Again the set of exceptions to the above 
limit is of the Lebesgue measure zero and it includes rational numbers, quadratic irrational etc. 



2 First experiment 

Let us define the sum of reciprocals of all Mersenne primes: 

B„ = £^> (8) 

n=l 

which can regarded as the analog of the Brun's constant, i.e. the sum of reciprocals of all twin 
primes: 




In 1919 Brun [5] has shown that this constant £>2 is finite, thus leaving the problem of infinity 
of twin primes not decided. Today's best numerical value is £>2 ~ 1.90216058, see [TB], [22] . 
Yet it is possible to prove that there is infinity of twins by showing that Brun's constant is 
irrational [27] (we believe it is even transcendental). 

Using PARI [21] we have calculated the sum of reciprocals of all known 47 Mersenne primes 
Bm with accuracy over 12 millions digits: 

B M = 0.5164541789407885653304873429715228588159685534154197 .... (10) 
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r- — 

Qk 



This number is not recognized by the Symbolic Inverse Calculator (http://pi4acim.uqam.ca) 
maintained by Simone Plouffe. The bar over Bm denotes the finite (at present consisting of 47 
terms) approximation to the full sum defined in ([8]). It is not known, whether there are Mersenne 
prime numbers with exponent 20996011 < p < 43112609 — currently confirmed by GIMPS is 
that 2 20996011 — 1 is the 40-th Mersenne prime A^4o — it is not known whether any undiscovered 
Mersenne primes exist between the 40th Ai^o and the 47th Mersenne prime .M47. We have taken 
12000035 digits of Bm — it means that we assume that there are no unknown Mersenne primes 
with p < 39863137. Using the incredibly fast procedure ContinuedFraction[-] implemented 
in Mathematica® we calculated the continued fraction expansion of Bm- The result was built 
from 11645012 partial denominators a\ = l,a 2 = l,o 3 = 14, . . . , 011645012 — 4. The n-th 
convergent P n (r)/Q n (r), see (|5f, approximate the value of r with accuracy at least 1/QkQk+i 
[T3l Theorem 9, p. 9]: 

111 

< 7T7r~ < W^ < cfi- 

From this it follows that if r is known with precision of d decimal digits we can continue with 
calculation of quotients a n up to such n that the denominator of the n-th convergent Q 2 n < 10 d . 
We have checked that Qn 64 5oi2 = 4.291385 x io 6000016 . 

The largest denominator was 09965536 = 716699617. We have checked correctness of the 
continued fraction expansion of Bm by calculating backwards from [0; 1, 1, 14, . . . , 4] the partial 
convergent. The Mathematica® has build in the procedure FromContinuedFraction [•] , but 
we have used our own procedure written in PARI and implementing the recurrence: 

n>l (12) 

with initial values 

P = o , Q = 1, Pi = o oi + 1, Qi = 01. (13) 

We have obtained the ratio of two mutually prime 6000018 decimal digits long integers (it 
means denominator was of the order io 6000018 and hence its square was smaller than io 12000035 ^ 
see eq.(JlT])): 

6000018 digits 



2216304109121123313251143869 . . . 2210 
4291385759849224534616716035 . . . 2813 

whose ratio had 12000033 digits the same as Bm- The decimal expansion of Bm — P/Q is of 
course periodic (recurring), see [9j Th. 135], but the length of the period is much larger than 
1.2 x 10 12 . According to the Theorem 135 from [9] the period r of the decimal expansion of 
Bm is equal to the order of 10 mod Q, i.e. it is the smallest positive r for which 

10 r = l(modQ). (14) 

Because Q being the product of all 47 Mersenne primes is of the order 3.509 . . . x io 86789810 ; we 
expect that the value of r is much larger than 10 12 . Another argument is that we got over 11 
500 000 partial quotients of the continued fraction of Bm — the numbers with periodic decimal 
expansions have only finite number of partial quotients different from zero. 
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From the sequence of partial quotients ai 
calculated running geometrical means 

K(n) = 



l,a 2 = I, a 3 = 14, 




Oii6450i2 = 4 we have 



(15) 



for n = 11, . . . , 11645012. The obtained numbers K{n) quickly tend to the Khinchine constant 
thus in Fig. [2] we have plotted the differences \K{n) — K\. The power fit to the values for 

~ n~ - 79 and it suggests that 
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1000 . . . 11645012 gives the decrease of the form \K(n) — K\ 



indeed lim^oo K{n) = K and thus Bm is irrational. Indices n for which the geometric means 
K{n) produce progressively better approximations to Khinchin's constant are: 



1, 3, 2, 16, 17, 21, 24, 26, 29, 412, 788, 1045, 369625, 369636, . . . , 5137093, 10389989; 



(16) 



the smallest value of \K(n) — K\ was 4.455957. . . x 10~ n . This sequence can be regarded as 
the counterpart to the A048613 at OEIS.org. 

Next we calculated running (i.e. for n = 11, 11645012) partial quotients P n /Q n and 
then the quantities L(n) = \/Q n , which for almost all irrational numbers should tend to the 



Khinchine-Levy constant. The behaviour of \[Qn is shown in Fig|3j Again we see that these 
quantities tend to the limit L; the fitting of the power-like dependence for n > 10 gives that 
\L(n) — L\ ~ 175.3977 -0 ' 92 . The shape of the plot in this figure is similar to the plot of \K{n)—K\ 
in Fig. [2} 

Both differences K(n) — K and L(n) — L have a lot sign changes for n < 11645012. Figures 
[4] and [5] present the plots of these differences together with the number of sign changes. 

The data presented in Figures [2] and [3] provide the hints that Bm is irrational and hence 
that there is infinity of Mersenne primes. But we are convinced Bm is in fact transcendental. 
In favor of this claim we recall here the result of A. J. van der Poorten and J. Shallit [25] that 
the following sum 

1111 1 

where F n are Fibonacci numbers, is transcendental. It is well known that the Liouville number 

1111 1 

is transcendental see [91 Theorem 192]. In Bm, assuming the Wagstaff conjecture, unfortunately 



the terms decrease slower: n! > 2 n > 2 e 



for n > 4 but faster than F„ 



V5 



where 



1, 2, . . . , n which are equal 



y = l±fl ~ 1.6180339887.... 

Let ip n (m) denotes the number of partial quotients a^ with k 
to m: 

' l Pn{ m ) — tt{^ '■ k < n and a^ = m}. 
Then the Gauss-Kuzmin theorem (for excellent exposition see e.g. [10, §14.3]) asserts that 

Mm) lo S ( X + 



lim 

n— >oo 



n 



log(2) 



(19) 
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for continued fractions of almost all real numbers. In other words, the probability to find the 
partial quotient = m is equal to log 2 (l + l/m(m + 2)). In Fig. [6] we present the plot of 
the ^Yi^o'is"' 1 f° r ^ ne continued fraction of Bm and m — 1, 2, . . . 1000 together with prediction 
given by the Gauss-Kuzmin theorem finding excellent agreement. 

Finally let us notice, that the number Bm computed with 12000035 digits is normal in the 
base 10, see Table I. 



TABLE I 

Illustration of the normality of Bm' the numbers in second row gives the number of digits 
0, 1, ... 9 appearing in the decimal expansion of Bm and the third row contains the ratio of 

numbers in second row divided by 12000035. 






1 


2 


3 


4 


1200553 


1199322 


1199420 


1200548 


1199397 


0.1000458 


0.0999432 


0.0999514 


0.1000454 


0.0999495 


5 


6 


7 


8 


9 


1198596 


1200876 


1200056 


1201757 


1199510 


0.0998827 


0.1000727 


0.1000044 


0.1001461 


0.0999589 



For comparison we have repeated the above procedure for artificial set of 47 numbers of 
the size corresponding to known Mersenne primes. We have simply skipped -1 in the Mersenne 
primes and using PARI we have computed with over 120000000 digits the sum: 

1 1 1 1 

2^ 2^ ^ ' ' ' ^ 2 42643801 ^ 2 43112609 

This number S is the ratio of the form 7 4 / /2 43112609 , where gcd(A, 243112609) = 1. From 
[HI §9.2] we know that S has terminating decimal expansion consisting of 43112609 decimal 
digits, thus calculating 12000000 digits of this sum makes sens as it does not contain recurring 
periodic patterns of digits. We have developed S into the continued fraction, what resulted in 
10550114 partial quotients. The calculated quantities for this case we denote with the subscript 
2: Q^in), K 2 (n), L 2 {n) to distinguish them from earlier experiment for Bm- We have calculated 
the running geometrical averages of the partial quotients K 2 (n) and the results are presented in 
Figure [7| Next we calculated 10550114 partial convergents P2(n)/Q 2 (n),n = 1,2, ... , 10550114 
and from them the quantities L 2 {n) = (Qzfa)) 5 which should tend to the Levy constant L. 
In Figure [8] the differences \L 2 {n) — L\ are plotted. Obtained plots are completely different 
from those seen in Figures [2] and [3] and they suggest K 2 {n) as well as L 2 {n) do not have the 
limit. In this artificial case we have encountered the phenomenon of extremely large partial 
denominators: there were a n of the order 10 70548 , 10 97732 and io 279910 . These large partial 
denominators are responsible for the smaller number of than for Bm, see (11). 



3 Second experiment 

Let us define the supposedly infinite and convergent continued fraction u_m by taking a n = M. n \ 



u M = [0;3, 7, 31, 127, 8191, 131071, 524287, 2147483647, . . .] (20) 
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Using all 47 Mersenne primes 3 , 7, 31, ... , 2 43112609 — 1 in a couple of minutes we have cal- 
culated um with the precision of 10000000 digits; first 50 digits of um ar e: 

u M = 0.31824815840584486942596202748140694243806236564 . . . (21) 

This number is not recognized by the Symbolic Inverse Calculator (http://pi.lacim.uqam.ca). 



Because l/Q%{u M ) « 1.84313 x icr 173579621 it follows from ([nj that theoretically it is possible 
to obtain the value of um from presently known 47 Mersenne primes with over 170,000,000 
decimal digits of accuracy. Of course um is the exception to the Khinchin and Levy Theorems 
in view of the very fast growth of u_m — see the Wagstaff [26] conjecture 0. 

There is a vast literature concerning the transcendentality of continued fractions. For ex- 
ample the continued fraction 

[0;2 1! ,2 2! ,2 3! ,...,2 n! ,...] (22) 

is transcendental, see Theorem 192], [231 Example 4]. 

The Theorem of H. Davenport and K.F. Roth [7] states, that if the denominators Q n of 
convergents of the continued fraction r = [ao; Oi, ci2, ■ ■ ■] fulfill 

y/\og(n) log(log(Q„(r))) 

hmsup — = oo (23) 

n n 

then r is transcendental. This theorem requires for the transcendence of r very fast increase of 



denominators of the convergents: at least doubly exponential growth is required for (23). The 
set of continued fractions which can satisfy the Theorem of H. Davenport and K.F. Roth is of 
measure zero, as it follows from the Theorem 31 from the Khinchin's book [13], which asserts 
there exists an absolute constant B such that for almost all real numbers r and sufficiently 
large n the denominators of its continued fractions satisfy: 

Qn(r) < e Bn . (24) 

The paper of A. Baker [4] from 1962 contains a few theorems on the transcendentality of 
Maillet type continued fractions pj)], i. e. continued fractions with bounded partial quotients 
which have transcencendental values. Besides Maillet continued fractions there are some spe- 
cific families of other continued fractions of which it is known that they are transcendental. 
In the papers [T7], [2] it was proved that the Thue-Morse continued fractions with bounded 
partial quotients are transcendental. Quite recently there appeared the preprint [6] where the 
transcendence of the Rosen continued fractions was established. For more examples see [3]. 



In the paper [lj B. Adamczewski and Y. Bugeaud, among others, have improved (23) to the 
form: 

Tf r log(log(Q n (r))) 

If hmsup^— — -^r, — j — r^\ =oc 25 
n n*' 6 Yogynp* log(log(n)J 

then r is transcendental. 

Assuming the Wagstaff conjecture Ai n ~ 2 2 " e mentioned in the Introduction we obtain 
that for large n 

Q n > 2 c2(n+1)e " 7 , c = - _ 7 X - = 2.101893933 . . . (26) 



s 



and thus the transcendence of um will follow from the Davenport-Roth Theorem (23): 

\A>g(n) log(log(Q n (r))) 



7? 



n — >■ oo. 



(27) 



We illustrate the inequality (26) in Figure |9| — the values of labels on the y-axis give an idea of 
the order of the fast grow of Q n (uM)' the largest for n = 47 is of the order Q47 = e 1 
2.32928 ... x io 86789810 , see also Table II. 



,1.9984.. .xlO 8 



TABLE II 

A sample of values of inverses of the squares of the n-th convergents of um giving an idea of 
the speed of convergence of [0; Aii, M2, ■ ■ ■ , -M-n] to um- 



Tl 


1 in 2 


Q 
O 


9 1 Q1 1 7Q7/1Q v 1 C\ — 6 


A 

4 


i.ozuoozoiy X 1U 


e 



i.yoo4ioyoy x ru 


P. 



i.i40(ooyoo x ru 


1 


4.lDoo04£)D0 X 1U 


Q 
O 


y.uoooyyo4z x ru 


9 


1.699990496 x 10" 95 






17 


9.32543401 x 10" 4439 


18 


1.38891910 x l(T 6375 


19 


3.81534516 x 10" 8936 


20 


4.67942175 x 10" 11599 






40 


4.50116310 x i(r 31553835 


41 


5.02100786 x 1 - 46025300 


42 


3.36434042 x 1(T 61657758 


43 


3.38166968 x io~ 79961861 


44 


2.17906011 x 10" 99578575 


45 


5.32688381 x i(T m949118 


46 


1.84595823 x lo^ 147623244 


47 


1.84313029 x 10" 173579621 



One of the transcendence criterion is the Thue-Siegel-Roth Theorem, which we recall here 
in the following form: 

Thue-Siegel-Roth Theorem: If there exist such e > that for infinitely many fractions 
A n /B n the inequality 



A„ 



< 



£24 

n 



n 



1,2,3,..., 



(28) 



holds, then r is transcendental. 



Let us stress, that e here does not depend on n — it has to be the same for all fractions 



A n /B n . We can test the criterion (28) for u_m using as the rational approximations A n /B n the 



convergents P n /Q n of the continued fraction (20). 
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In [23] J. Sondow has given the estimation for e appearing in r.h.s.of (28); namely he proved 
that for irrational numbers with continued fraction expansion [ao',a>i,a2, ■ ■ •] and convergents 

Pn/ Qn- 

loga n+1 

e < hmsup - — — — . (29) 

n— >oo lOg (%n 

Let us denote 5 = 2 + e. From the Wagstaff conjecture we obtain that the exponent S of B s 



appearing in on the r.h.s. of (28) should be of the order 



5^2 + 2 e ~ 7 - 1 = 2.47477 . . . (i.e. e « 0.47477 . . .) (30) 



implying transcendence of um- in Fig, 10 we present actual values of 5(«»;n) = — log \um — 
Pn/Qn\/l°g(Qn) for n = 3, 4, . . . , 45 and indeed the values oscillate around l+2 e 7 = 2.47477 
First we have calculated um using all 47 Mersenne primes with accuracy 140000000 digits and 
for n — 3, 4, . . . , 45 we have calculated convergents P n /Q n and next the differences \um~ P n /Qn\ 
with accuracy (see Table II), from which we determined the 5(uM\ n )- The arithmetic 



average of 43 values 5(w^;n) is 2.5002..., quite close to the estimated value (30). It took 



a few months of CPU time to collect data presented in Fig. 10 It took 12 days of CPU 
time on the AMD Opteron 2700 MHz processor to collect data for n < 40; the point n = 40 
needed precision of almost 40,000,000 digits, as \u M - P40/Q40I = 1-5033 x 1(T 38789567 , while 
l/Q\ = 4.501 ... x 10- 31553835 . To calculate the difference \u M - P n /Q n \ for n = 41, 42, 43 the 
precision of 100000000 digits was needed. For n = 44 and n = 45 the difference \um — P n /Q n \ 
was calculated with the precision 130000000 digits (see Table II for n = 44 and n = 45) and it 
took about one month of CPU time for each point. 
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Figure 1: The plot of log log (.A/In) arid the Wagstaff conjecture Qlh . The fit was made to all 
known M n and it is 0.3854n + 0.6691, while ne" 7 log(2) - loglog(2) « 0.3892n + 0.3665. The 
rather good coincidence of log log (.M n ) and Q is seeming, as to get original .M n 's the errors 
are amplified to huge values by double exponentiation. 




Figure 2: The plot showing the distance to K of the running geometrical averages K( 
(aia^ • • • a n ) 1//n for the continued fraction of Bm- 
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Figure 3: The plot showing the distance to L of the (Q(n)Y^ n obtained from the partial 
convergents of the continued fraction of Bm for to = 11, ... , 11645012. 
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Figure 4: The plot of K(n) in black approaching the Khinchine constant K = 2.685452 ... (in 
red) with values presented on left y-axis. In green are presented numbers of sign changes of 
K(n) — K up to n — the right y-axis is for this plot. 
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Figure 5: The plot of L(n) in black approaching the Khinchine-Levy constant L = 
3.27582291872 ... (in red) with values presented on left y-axis. In green are presented numbers 
of sign changes of L(n) — L up to n — the right y-axis is for this plot. 




Figure 6: The plot of the measured for the continued fraction of Bm probability to find the 
partial quotient = m for the continued fraction of Bm- 
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Figure 7: The plot showing the distance to K of the running geometrical averages K 2 (n) for 
the continued fraction of S for n — 11, ... , 10550114. 
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Figure 8: The plot showing the distance to L of the (Q2(n)Y^ n obtained from the partial 
convergents of the continued fraction of S for n = 11, . . . , 10550114. 
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Figure 9: Illustration of the inequality (26) for 3 < n < 47. Although the last points seem 
to coincide in fact Q 47 = 2.32928... x io 86789810 , while T 2 ^~~' = 1.21513... x io 82034318 - 
hundreds thousands orders of difference! 
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Figure 10: The plot of — log \u M — P n /Q n \/ \og(Q n ) fluctuating around the estimation (30). 



